Projecting the Kondo Effect: Theory of the Quantum Mirage 
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A microscopic theory is developed for the projection (quantum mirage) of the Kondo resonance 
from one focus of an elliptic quantum corral to the other focus. The quantum mirage is shown to 
be independent of the size and the shape of the ellipse, and experiences Af/4 oscillations (A_p is the 
surface-band Fermi wavelength) with an increasing semimajor axis length. We predict an oscillatory 
behavior of the mirage as a function of a weak magnetic field applied perpendicular to the sample. 

PACS numbers: 72.15.Qm, 61.16.Ch, 72.10.Fk 

In a recent experiment, Manoharan et al. [jlj used an 
elliptic quantum corral to project the image of a Kondo 
resonance over a distance of tens of angstroms, from one 
focus of the ellipse to the other focus. By placing a mag- 
netic Co atom at one focus of the ellipse and measuring 
the tunneling current to a close-by scanning tunneling 
microscope (STM) tip, a distinctive Kondo resonance was 
seen in the I-V curve when the tip was brought directly 
above the Co adatom. Remarkably, a similar Kondo sig- 
nature was observed when the tip was placed above the 
empty focus, indicating coherent refocusing of the spec- 
tral image by the surrounding corral. This should be 
contrasted with STM measurements of isolated magnetic 
adatoms on open surfaces [pH, where a limited spatial 
extent of ~ 10A was observed for the Kondo effect. 

Semiclassically, one can attribute this refocusing phe- 
nomena to the property that all classical paths leaving 
one focus of the ellipse bounce specularly off the perime- 
ter and converge onto the second focus with the same ac- 
quired phase Q| (see Fig. 1). However, this simple picture 
does not explain the quantitative features of the experi- 
ment. For example, the complex interference patterns in 
the dl/dV difference map throughout the ellipse, or the 
Xp/i oscillations of the mirage with an increasing semi- 
major axis length a (Xp is the Fermi wavelength). Ex- 
planation of these features requires a quantitative theory, 
which is the objective of the present Letter. 

Starting with a microscopic picture of Kondo scatter- 
ing off the Co adatoms we obtain good qualitative and 
quantitative agreement with the experiment. We estab- 
lish a remarkable feature of the quantum mirage which, 
aside from the A^/4 oscillations mentioned above, is inde- 
pendent of the size and the shape of the ellipse, provided 
the ellipses is not too small. In particular, there is no 
dependence on the ellipse eccentricity £ , see Fig. 1. In 
the presence of a weak perpendicular magnetic field, we 
predict an oscillatory behavior of the quantum mirage as 
a function of the magnetic flux encircled by the ellipse. 

The Cu(lll) surface has a band of surface states, 
which acts as a two-dimensional electron gas. The sur- 
face band starts 450meV below the Fermi energy, and 
has a Fermi wave number of hp 1 ~ 4.75A. When a Co 
adatom is placed on the surface, it scatters both the sur- 
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FIG. 1. An illustration of an elliptic quantum corral and 
the classical trajectories associated with the quantum mirage. 
The ellipse is characterized by a semimajor axis length a and 
an eccentricity £. The mean distance between adjacent atoms 
forming the ellipse perimeter is d. 



face electrons and the underlying bulk electrons. As re- 
cently shown by Ujsaghy et al. || for Co on Au(lll), one 
may model the Co adatom by an effective nondegener- 
ate Anderson impurity ||, characterized by an effective 
energy level ed, an on-site repulsion U, and two hybridiza- 
tion matrix elements t s and tb to the underlying surface 
and bulk conduction electrons. In this manner, each of 
the Co atoms forming the ellipse in the experiment of 
Manoharan et al. |l[] acts as an Anderson impurity, as 
does the adatom placed inside the ellipse. 

Denoting the creation of a surface-state and a bulk 
conduction electron by ct and at a , respectively (here 

k labels a two-dimensional surface vector while q is a 
three-dimensional vector), we model the system by the 
Hamiltonian Ti = 7i surf + 7i b uik + X«=o ^im P (-Ri)> where 

Wsurf = J2 e fe c L C £r and Hhum = 51 E 1 a l* a fc 



describe the free surface and bulk conduction bands, re- 
spectively, and 

T-Cim P {Ri) = e d ^2 d L d ^ + u d \i d n d ii d il 

a 

+ Y, { l AMRi) + Ud\ a Xo{Ri) + h.c.} 

(X 

describes a Co adatom at point R% on the surface. Here 
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d\ a creates a localized Co electron at site R4 (i = for 
the inner adatom and i = 1, • • • , N for the perimeter 
adatoms), while ip rT (Ri) and x a {Ri) annihilate, respec- 
tively, a surface and a bulk conduction electron at site Ri. 
For simplicity, we have taken the different adatoms to be 
identical, and neglected any momentum dependence of 
t s and tb- In what follows, we shall mainly be interested 
in the case where the inner adatom is located at the left 
focus, i.e., Rq = —R in the notations of Fig. 1. 

Consider now an STM tip placed directly above the 
surface point r. If the tip couples predominantly to the 
underlying surface-state electrons at f, then the differ- 
ential conductance for the current through the STM tip 
measures, up to thermal broadening, the local surface- 
electron density of states at point r, p(f,e). For an iso- 
lated Co impurity, p(r, e) depends on the Kondo scat- 
tering from the impurity as described, e.g., in Ref. 0. 
For the multiple-impurity configuration considered here 
there are two main modifications: (i) There are multiple 
scattering off the different Co adatoms; (ii) Intersite cor- 
relations alter the Kondo scattering off each Co adatom. 
Due to the relatively large distance between Co atoms 
(~ 10A for neighboring atoms on the ellipse perimeter), 
we expect the latter effect to be small, and therefore ne- 
glect it hereafter. 

Since we are mostly interested in the effect of the Co 
atom placed inside the ellipse, we distinguish it from 
the other Co atoms on the perimeter of the ellipse. 
Neglecting intersite correlations, p(r, e) takes the form 
p(f, e) = p(f, e) + Sp(f, e), where 



P(r, e) 
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lm{G(r,f;e)} 



(1) 



is the density of states of an empty ellipse (i.e., in the 
absence of the inner adatom), and 



6p(f, e) 



: Im{^G(f, R Q ; e)G d (e)G(R ,f; e)} (2) 



is the additional contribution due to the extra Co atom 
at -Ro- Here G(r, f '; e) is the retarded Green function of 
the surface electrons for an empty ellipse, and G<j( e ) is 
the fully dressed retarded Green function of the d elec- 
trons of the inner adatom. Note that in writing Eq. (||) 
we have assumed that the three-dimensional propagation 
of bulk electrons between different Co sites on the surface 
is small compared to the two-dimensional propagation of 
the surface electrons. This assumption is quite reason- 
able considering that the three-dimensional propagation 
near the surface decays as 1/r 2 , compared to l/y/r for 
the two-dimensional surface propagation . 

Experimentally, 5p(r, e) is extracted by first measuring 
the local density of states of the empty ellipse, and then 
subtracting it from the measured density of states with 
the extra Co atom. To compute Sp(f,e), one needs to 
evaluate G(r, r '; e), which is our next goal. To this end, 



we introduce the Nx N matrix = (l—Si > j)Go(Ri, Rj), 
along with the two vector quantities, Uj = Ga(r, Ri) 
and m — Go(Ri,r r ). Here Go(r,r') is the free surface 
Green function without the corral, and i and j run over 
1, • • • , N. Using these quantities, G(f, f; e) is compactly 
expressed as 



G(f,r';e) = G (r,f') + 



N 
i,3'=l 
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(3) 



where T(e) = t 2 Gd{e) is the surface-to-surface compo- 
nent of the conduction-electron scattering T-matrix at 
each Co site. Again, Eq. (^) omits the intersite correla- 
tions and the bulk propagation between different Co sites. 
Finally, for temperatures and energies below the Kondo 
temperature, the Kondo part of Gd(e) may be well ap- 
proximated [|| by the Lorentzian form Zk/(c— ep+iTx ), 
where Tk is the Kondo temperature, ep is the Fermi en- 
ergy, and 



Z h 



■np s t 2 s + npbt? 



(4) 



is the corresponding weight. Here p s and pi, are the sur- 
face and bulk density of states at the Fermi level. 

A key parameter that enters the quantum mirage is 
the ratio of scattering rates 



t 



■Kpstl 



--. (5) 

irpstj + np b tl 

Physically, t represents the probability that a surface- 
state electron impinging on a Co adatom will be scattered 
to a surface-state electron rather than a bulk electron. 
Hence f is a measure of the in-elasticity of the scattering 
of surface waves from the Co impurities. In the theory 
of Heller et al. || for the standing waves formed in a 
quantum corral, t is found to be 1/2. Hereafter we shall 
use the same value for t. 

In Fig. 2 we depict Sp(r, €f) for various configurations 
of the Co atoms, as measured by Manoharan et al. [0. 
The upper panels correspond to an ellipse with eccentric- 
ity E = 0.786 and 34 adatoms (ellipse b of Ref. 0), while 
the lower panels correspond to an ellipse with eccentricity 
£ = 0.5 and 36 adatoms (ellipse a in Ref. M). The quan- 
tum mirage is clearly seen in each of the left two panels, 
where an additional adatom has been placed at the left 
focus of the ellipse. For both ellipses, there is a strong 
signal in the tunneling density of states right above the 
right focus, in accordance with the experimental data. 
By contrast, the quantum mirage disappears when the 
additional adatom is places off the focus, as shown in the 
right two panels. These results are in good agreement 
with the experimental measurements, reproducing even 
fine details of the experimental patterns. 

The results of Fig. 2 were obtained from Eqs. (||) and 
(^), by setting t = 1/2 and approximating Go with the 
free two-dimensional Green function: 
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FIG. 2. Contour plot of 5p(f, ep) at the Fermi level assum- 
ing t — 1/2. The upper panels correspond to an ellipse with 
eccentricity £ — 0.786 and 34 adatoms, while the lower panels 
to an ellipse with £ = 0.5 and 36 adatoms. For both ellipses, 
a strong signal appears at the right focus of the ellipse when 
an additional adatom has been placed at the left focus. This 
quantum mirage disappears when the additional adatom is 
places off the focus, as shown in the right two panels. 



G (f,r") = -iK Ps H£\k\r-r'\). 



(6) 



Here k is the wave number, and (x) is the Han- 
kcl function of zeroth order, which for x 3> 1 takes the 
asymptotic form 



2 (. ,7T\ 

exp yix — i—j 



(7) 



The good agreement between our calculations with 
t = 1/2 and the experimental data indicates that the 
number of scattering events that a particle undergoes be- 
fore leaving the surface bounded by the ellipse is small. 
This suggests the possibility of calculating the quantum 
mirage at the right focus at the Fermi level, Sp(R,ep), 
using perturbation theory in t. Thus, to linear order in 
t, the Green function between the left and the right foci 
is given by 



G( — R, R] ep) — Go ( — R, R) + G\ ( — R, R) + • • • 
where 



Gq{—R, R) 



-ipi 



ilSh-pa—i^ 



Ekpa 



(8) 



(9) 



is the contribution of the direct path connecting the two 
foci (illustrated by the dashed line in Fig. 1), and 



G X {-R,R) = — ^Go(-R,Rj)G (Rj,R) (10) 



n 



iirp. 



comes form all trajectories which scatter from a single Co 
adatom on the ellipse perimeter (solid lines in Fig. 1). 

Conventionally, the second term G\ is smaller than Go 
for several reasons. First, each scattering event is inelas- 
tic, and therefore introduces a reduction factor of t. Sec- 
ond, the scattered orbits are longer. Third, the various 
orbits have generally different lengths, and therefore their 
corresponding phases add up incoherently. The situation 
is quite different in the present case. Due to the defining 
property of the ellipse, the length of all scattered orbits 
between the two foci is precisely the same, and hence 
their contributions add up coherently. Consequently, the 
second term in Eq. (||) takes the form 



TV 



2t 



i2kp-a—ii 



(11) 



where 2a is the length of each orbit, and £ij and 
are the distances between the impurity at Rj and the 
right and left foci, respectively (see Fig. 1). Finally, we 
approximate the sum in Eq. (|llj) by an integral: 



N 

E 



2 J 



d f ds VUs)W) 



(12) 



where s denotes the coordinate along the ellipse contour, 
and d is the mean distance between adjacent adatoms. 
The result of the integral is independent of the eccentric- 
ity of the ellipse, and is simply 27r. Thus, the contribution 
of orbits scattered from a single perimeter adatom is 



Gl (-R,R)^ Ps ^ 2k - a . 

KfCL 



(13) 



Comparing Go{~R,R) and G\{—R,R) at the Fermi 
energy, one sees that the leading contribution to the 
quantum mirage comes from Gi 0, provided 



d 167rf 2 
£a kpd 



(14) 



Substituting the experimental parameters, d ~ 10A, 
a ~ 70A, and k F ~ 1/4. 75A" 1 , and setting t = 1/2, 
it is straightforward to verify that Eq. (|lj) holds for all 
ellipses with eccentricity 0.05 < £ < 1. Furthermore, 
Eq. (|l4|) is always satisfied for sufficiently large ellipses, 
provided the mean distance between adjacent adatoms is 
kept fixed (a^> d). 

Neglecting the contribution of the direct path, 
Go(—R,R), and using Eqs. (0) and ( |l3| ) with r = R, 
the resulting local density of states at the Fermi energy 
takes the form 



5p(R,€ F ) ~ p s 



16i J 



■ cos(4fcp-a). 



(15) 



' (k F d) 2 

The main feature of the above result is the robustness of 



the quantum mirage: As long as condition (14) is satis- 
fied, the amplitude of the mirage is independent of the 
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size of the ellipse, a, and its eccentricity, £. Rather, the 
amplitude is determined by t, which characterizes the in- 
elasticity of the scattering of surface waves from adatoms, 
and the dimensionless mean distance between adjacent 
adatoms along the ellipse, kpd. The oscillations of the 
mirage as function of a are indeed periodic with a period 
of Ajr/4, as seen experimentally Q. 

Next we consider the effect of a weak uniform mag- 
netic field, B, applied perpendicular to the surface. As 
shown below, the quantum mirage experiences a distinc- 
tive oscillatory behavior as a function of the magnetic 
field, which depends on the size and the shape of the el- 
lipse. Here we assume that the ellipse is sufficiently large 
and that the magnetic field is sufficiently weak so that 
(i) Zeeman splitting of the Kondo resonance can be ne- 
glected, and (ii) the cyclotron radius of the conduction 
electrons is much larger than the ellipse size. 

Under these circumstances, the main effect of the mag- 
netic field is to introduce an additional Aharonon-Bohm 
phase to the contribution of each path. This phase is 
the magnetic flux encircled by the orbit, measured in 
units of the quantum flux, </>o = hc/e. Here h is Planck's 
constant, c is the velocity of light, and e is the elec- 
tron charge. To compute Gi(—R,R) we fix the gauge 
by calculating the flux encircled by the path which goes 
from R to Rj, to —R, and then back to R along the 
scmimajor axis of the ellipse. Accordingly, the sum of 
Eq. (|l|) is modified to Ej(Cij6j)~ 1/2 e l27r¥,j/0 °, where 
ipj is the flux of each trajectory. Using the continuum ap- 
proximation of Eq. ( |i~2] ) this sum gives 2ttJ q (2£ AB/(j) ), 
where Jo(x) is the Bessel function of zeroth order, and 
A = ita 2 \/l — £ 2 is the area of the ellipse. The quantum 
mirage is, thus, modified according to 

5 P (B) ~ 5p(0) J 2 (2£^) , (16) 

where 5p(0) is the zero- field result of Eq. (|l5|). Note that, 
for a given a, the sensitivity to a magnetic field is largest 
for £ = 1/V2. 

Finally, we discuss the role of imperfections in the el- 
lipse. Clearly, a combined effect of many coherent tra- 
jectories is very sensitive to imperfections and dephas- 
ing. At 4K, dephasing effects due to electron-electron 
and electron-phonon interactions are negligible over dis- 
tances of the order of hundreds of angstroms. In what 
follows we show that effect of imperfections is small too. 

The main source of imperfections in the ellipse comes 
from the position of the adatoms forming the ellipse. 
These are constrained to sit on a triangular lattice im- 
posed by the underlying Cu(lll) surface. Consequently, 
the lengths of the orbits contributing to Sp are not those 
of an ideal ellipse. To estimate the effect of the devia- 
tions, we consider a random distribution of the trajectory 
lengths, and average over the distribution. We first notice 
that each contribution to 5p is composed of 4 segments: 



£ij £2, £3, and £4, as illustrated in Fig. 1. Each of these 
segments is regarded as an independent random vari- 
able uniformly distributed in the range (£ — 6/2, £ — 6/2), 
where £ is the exact distance from the focus to the ellipse 
boundary, and 6 is the triangular lattice spacing. The to- 
tal length 77 = J2 n =i 4 ^ s > therefore, approximately a 
Gaussian random variable with mean 4a and variance 
6 2 /3. Averaging the cosine term in Eq. (p"5|), (cos(kpri)}, 
the disorder-averaged value of the quantum mirage is re- 
duced by a factor of approximately Q — exp(— k F b 2 /6). 
Substituting the experimental values kp — 1/4.75A _1 
and 6 = 2.55A, one finds Q « 0.95, meaning that the 
effect of imperfections is negligible. 

In conclusion, in this Letter we have studied the phe- 
nomenon of quantum mirage, and clarified its relation to 
the classical orbits of a particle in an ellipse. Our the- 
ory also predicts a distinctive behavior of the quantum 
mirage in the presence of a perpendicular magnetic field, 
which could be tested experimentally. Finally our ap- 
proach clearly shows that the phenomenon of quantum 
mirage is not unique for magnetic adatoms. It will also 
appear for nonmagnetic atoms (on the ellipse perimeter 
or at its focus) with strong scattering, e.g. adatoms with 
resonant tunneling states at the Fermi level. 
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